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CALCULATION OF FLUTTER CHARACTERISTICS FOR FINITE -SPAN 
SWEPT OR UNSWEPT WINGS AT SUBSONIC AND SUPERSONIC 
SPEEDS BY A MODIFIED STRIP ANALYSIS 
By E. Carson Yates, Jr. 

SIMMARY 


A method has been developed for calculating flutter characteristics 
of finite-span swept or unswept wings at subsonic and supersonic speeds. 
The method is basically a Rayleigh type analysis and is illustrated with 
uncoupled vibration modes although coupled modes can be used. The aero- 
dynamic loadings are based on distributions of section lift-curve slope 
and local aerodynamic center calculated from three-dimensional steady- 
flow theory. These distributions are used in conjunction with the 
"effective" angle-of-attack distribution resulting from each of the 
assumed vibration modes in order to obtain values of section lift and 
pitching moment. Circulation functions modified on the basis of loadings 
for two-dimensional airfoils oscillating in a compressible flow are 
employed to account for the effects of oscillatory motion on the magni- 
tudes and phase angles of the lift and moment vectors. 

Flutter characteristics have been calculated by this method for 
12 wings of varying sweep angle, aspect ratio, taper ratio, and center- 
of-gravity position at Mach numbers from 0 to as high as 1.75- Compari- 
sons of the results with experimental flutter data indicate that this 
method gives generally good flutter results for a broad range of wings . 


INTRODUCTION 


Much of the difficulty encountered in attempting to predict flutter 
characteristics for finite-span swept and unswept wings at subsonic and 
supersonic speeds results from inadequate representation of the distri- 
butions of oscillating aerodynamic loads on such wings. For both sub- 
sonic and supersonic speeds a number of methods exist for evaluating 
three-dimensional oscillating loads (refs. 1 to 21, for example). These 
methods involve varying degrees of rigor, but all are characterized by 
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the relatively extensive amount of computation required. In all of these 
procedures it is necessary to recalculate the loading with each change of 
reduced frequency. This fact further increases the amount of computation 
required because in flutter prediction the reduced frequency at flutter 
is not usually found directly. Because the calculations are complex and 
lengthy and because many of the procedures have not been proved in gen- 
eral application, the use of these methods in flutter prediction has been 
limited. 

A procedure commonly used in the solution of practical flutter prob- 
lems involving finite wings is a modal-type analysis similar to that 
employed by Barmby, Cunningham, and Garrick for swept wings (ref. 22) 
and by Smilg and Wasserman for unswept wings (ref, 23). These methods, 
as presented in references 22 and 23 , employ two-dimensional incompres- 
sible aerodynamic forces and moments and thus do not take Into account 
the aerodynamic effects of finite span and compressibility. 

The present report presents an approximate method of flutter cal- 
culation based on a simplified representation of the three-dimensional 
aerodynamic loading which is shown to be applicable to a wide variety 
of wing plan forms at both subsonic and supersonic speeds. The present 
method is also based on a modal analysis, but the aerodynamic effects 
of finite span, taper, and compressibility are accounted for by utilizing 
modified aerodynamic loadings based on spanwise distributions of section 
lift-curve slope and local aerodynamic center calculated from well-known 
subsonic (ref. 24) or supersonic (refs. 25 and 26 ) three-dimensional 
steady-flow theory for flat, rigid wings. The distributions of section 
lift and pitching moment on oscillating flexible wings are obtained by 
employing these distributions of lift-curve slope and aerodynamic center 
for flat rigid wings in conjunction with the "effective" angle-of-attack 
distribution resulting from oscillation of the wing in each of the assumed 
vibration modes. The effect of oscillatory motion on the magnitudes and 
phase angles of the lift and moment vectors is represented approximately 
by modifying the familiar circulation functions of Theodorsen by utilizing 
aerodynamic flutter coefficients given by Jordan (ref. 27) for two- 
dimensional airfoils oscillating in subsonic or supersonic flow. A 
detailed description of the procedure for making flutter calculations is 
given in the appendixes. 

By representing the oscillating aerodynamic loads in this manner 
the necessity of recalculating the load distributions for each value of 
reduced frequency is avoided, since only the modified circulation func- 
tions vary with frequency, and these in turn are assumed not to vary 
along the span. The bending and twisting deformation of individual wing 
sections is taken into account only in terms of the "effective" angle of 
attack and is assumed not to affect distributions of lift-curve slope 
and aerodynamic center. This procedure is equivalent to neglecting the 
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influence of deformation on the lift -producing capacity of a given wing 
section. 

Flutter characteristics have been calculated by the method developed 
herein (using three vibration modes) for wings with sweep angles from 0° 
to 52 . 5 °, aspect ratios from 2.4 to 7.4, taper ratios of 0.6 and 1.0, and 
center-of -gravity positions between 34 percent chord and 59 percent chord. 
The results are compared herein with experimental data obtained in the 
Langley 26 - inch transonic blowdown tunnel (refs. 28 to 31) sJid in "the 
Langley 9- by 12-inch supersonic blowdown tunnel (ref. 32). 


SYMBOLS 


A 


a 


ac 


ac n 


b 

b r 

s 


B 



aspect ratio of full wing including fuselage intercept 

aspect ratio of wing considering side of fuselage as a reflec- 
tion plane (twice the panel aspect ratio) 

nondimensional distance from midchord to elastic axis measured 
perpendicular to elastic axis, positive rearward, fraction 
of semichord b 

nondimensional distance from leading edge to local aerodynamic 
center (for steady flow) measured streamwise, fraction of 
streamwise chord, a 

nondimensional distance from midchord to local aerodynamic 
center (for steady flow) measured perpendicular to elastic 
axis, positive rearward, fraction of semichord b 

semichord of wing measured perpendicular to elastic axis 

semichord of wing measured perpendicular to elastic axis at 
spanwise reference station tj = 0.75 

span of wing panel considering side of fuselage as a reflec- 
tion plane 

ratio of local semichord b to reference semichord b r meas- 
ured perpendicular to elastic axis, b/b j. 

complex circulation function, F + iG 

local lift-curve slope for a streamwise section in steady flow 
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C l 


a,n 


c ma 


C p 

F 

f h 

*0 

G 

g 

h 

la 

i 

knr 


l 

M 

Ma 


m 


P 

Q 


local lift-curve slope for a section perpendicular to the 
elastic axis in steady flow 

derivative with respect to angle of attack of local pitching- 
moment coefficient measured about the leading edge of a 
streamwise section 

local lifting-pressure coefficient 


circulation function which modifies in-phase load components 

deflection function of wing in bending mode 

deflection function of wing in torsion mode 

circulation function which introduces out-of -phase load 
components 

structural damping coefficient for wing (Subscript a denotes 
torsional mode; subscript h denotes bending mode.) 

local vertical translational displacement of wing at elastic 
axis 

mass moment of inertia of unit length of wing about e las tic 
axis 

i FT 

reduced frequency based on the spanwise reference station 
(l = 0-75) and on velocity component normal to elastic axis, 
b r a>/v n 


length of exposed wing panel measured along elastic axis 
Mach number 

oscillatory moment about elastic axis per unit length of wing, 
positive leading edge up 

mass of wing per unit length measured along elastic axis 

oscillatory lift per unit length of wing along the elastic 
axis, positive downward 

downwash expression defined by equation ( 5 b) 


H 
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r a 

t 

V 

Vr 

v 

X 

x' 

t 

y 

a 

P 

K 

A 

A 

X P 

1 

0 

P 


nondimens ional radius of gyration of wing about elastic axis. 



time 

flutter speed, measured parallel to free stream (experimental 
values or values calculated by the method of this report) 

calculated reference flutter speed obtained by using CS = 2 jc 

XI 

and ac n - - — 
n 2 

free-stream velocity 

streamwise coordinate measured from leading edge of wing root 

nondimens ional coordinate from midchord measured perpendicular 
to elastic axis, positive rearward, fraction of semichord b 

nondimensional distance from elastic axis to local center of 
gravity measured perpendicular to elastic axis, positive 
rearward, fraction of semichord b 

distance along elastic axis measured from wing root, lr\ 

angle of attack 

for M > 1; yi - M 2 for M < 1 
wing section mass-density ratio, jtpb 2 /m 
sweep angle; positive for sweepback 

taper ratio of full wing including fuselage intercept 
taper ratio of exposed wing panel 

nondimensional coordinate (either spanwise or along elastic 
axis) measured from wing root, fraction of exposed panel 
span s or fraction of wing length l 

local torsional displacement of wing measured about elastic 
axis 

air density 
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cr local bending slope of elastic axis, dh/dy' 

t local rate of change of twist, d0/dy r 

0) circular frequency of vibration 

o*x circular frequency of first uncoupled torsional vibration mode 

of wing measured about elastic axis 

cjq^ circular frequency of uncoupled bending vibration mode of wing 

(subscripts 1 and 2 denote first and second bending modes) 

£ nondimensional streamwise coordinate measured from leading edge 

of wing root, fraction of exposed panel span s 

Subscripts: 

c /4 quantities associated with the wing quarter-chord 

ea quantities associated with the wing elastic axis 

C circulation functions obtained from the oscillatory aerodynamic 

coefficients given in reference 27 for two-dimensional com- 
pressible flow 

LE quantities associated with the wing leading edge 

M quantities associated with the Mach lines originating from wing 

root or tip 

n quantities associated with wing sections normal to the elastic 

axis 

I circulation functions obtained by Theodorsen in reference 33 

for two-dimensional incompressible flow 

TE quantities associated with the wing trailing edge 

Dots over symbols denote derivatives with respect to time. 
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DESCRIPTION OF THE METHOD 


General 

The procedure for flutter calculation used in this report is bas- 
ically a Rayleigh, or modal-type, analysis and is illustrated herein with 
uncoupled vibration modes although coupled modes can be used. (The use 
of uncoupled modes in flutter calculations is discussed in detail in 
refs. 22 and Jk.) The flutter modes of the wings studied in this inves- 
tigation are represented by the first and second bending and the first 
torsional vibration modes of uniform cantilever beams. All deformations 
are considered to be made up of vertical bending of an approximately 
straight elastic axis and rotation about that axis. The wing root is 
treated as though it were clamped along a line normal to the elastic axis 
and passing through the intersection of the elastic axis and the root 
chord. The dynamical equations involved in this type of analysis are 
obtained from Lagrange's equations of motion in which the vibration modes 
are used as generalized coordinates. These dynamical equations repre- 
senting the balance between elastic, inertial, and aerodynamic loads are 
derived in appendix A and are obtained (for the simple case of one bending 
mode and one torsion mode) in the form 


^ * *> - 


X 


o' -(f) 


h + 


-bp 


x: m 


*a f h f 8 d y’ 


9 - 


*pb r 2 a£ J 0 


f l Pf h dy 1 = 0 
Jo 


( 1 ) 


and 


r 1 



~ w f 0 K^) Wedy ' 

h + ■ 

t^ 2 

where h and 0 

are 

as 


5-1 1 * «-) ' 


X 




Kpbp 2 ^ 


r 

Jo 


Maf 0 dy' = 0 


(2) 


i) . These same 
The values of 


equations in a different form were used in reference 22 
all geometrical, structural, and aerodynamic quantities to be used in 
these equations are those values associated with sections normal to the 
elastic axis. 


The innovations of the present method consist of alterations in the 
expressions for section lift P, pitching moment Ma, and complex circu- 
lation function C = F + iG in order to approximate the aerodynamic 
effects of finite span, taper, and compressibility. The section lift P 
and pitching moment M<x are expressed in terms of arbitrary section lift- 
curve slope and aerodynamic center which are assumed to vary along the 
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span of the wing. For any particular value of free-stream Mach number , 
the spanwise distributions of lift-curve slope and aerodynamic center 
are calculated from well-known steady- state aerodynamic theory for flat 
rigid wings. The spanwise distributions of the lift and moment on the 
deforming wing are then found by using the aforementioned values of 
static section lift-curve slope and aerodynamic center in conjunction 
with the '’effective” angle-of-attack distribution resulting from oscil- 
lation of the wing in each of the assumed vibration modes. 1 The values 
of lift and moment thus obtained account approximately for finite span, 
taper, compressibility, and deformation shape of the wing. However, it 
is also necessary to take into account the effect of oscillatory motion 
on the magnitudes and phase angles of the lift and moment vectors. In 
the present method this is done approximately by utilizing circulation 
functions (analogous to the familiar F and G functions of Theodorsen 
(refs. 35 and 55)) which are modified on the basis of aerodynamic flutter 
coefficients given by Jordan (ref. Zj) for two-dimensional airfoils oscil- 
lating in subsonic or supersonic flow. In the application of the circu- 
lation functions thus obtained, the Mach number normal to the leading 
edge is employed. 

Formulating the aerodynamic forces and moments in this manner implies 
the following assumptions: 

(1) The bending and twisting deformation of individual wing sections 
Is accounted for in terms of the "effective” angle of attack only. The 
effect of relative deformation on section lift -curve slope and aerodynamic 
center can be neglected. Camber deformation of sections normal to the 
elastic axis is not considered. 

(2) The effect of oscillatory motion on the magnitude and phase 
angles of the section lift and moment vectors is the same for each wing 
section and may be represented by modified circulation functions associ- 
ated with the Mach number component normal to the leading edge. 

In view of the use of static lift-curve slopes and aerodynamic cen- 
ters, application of this method at high values of reduced frequency 
would be open to question. At low to moderate reduced frequencies, how- 
ever, the approximation should be reasonable. 

In the remaining sections of this description of the method are 
discussed the alteration of section lift P and pitching moment Mq, 
by the introduction of static three-dimensional section lift-curve slopes 


^The "effective” angle of attack is the downwash resulting from the 
motion divided by the component of free-stream velocity normal to the 
elastic axis. 
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and. aerodynamic centers, the calculation of these static aerodynamic 
parameters, and the evaluation of the complex circulation function C 
by utilizing two-dimensional subsonic or supersonic oscillating-airfoil 
theory . 

A detailed description of the flutter calculation procedure is given 
in appendix B, and expressions for the elements of the final flutter 
determinant are given in appendix A. 


Expressions for Section Lift and Pitching Moment 

In formulating the expressions for section lift and pitching moment 
the following basic assumption is made: The flow over wing sections nor- 

mal to the elastic axis consists of a quasi-two-dimensional noncirculatory 
flow plus a circulatory flow in which the circulation is fixed by the 
component of free-stream velocity normal to the elastic axis in conjunc- 
tion with downwash distributions along chord lines normal to the elastic 
axis (rather than by the free-stream velocity and downwash distributions 
along streamwise chord lines). In contrast to the method of reference 22 
the present method does not consider the circulatory flow to be two- 
dimensional and incompressible in nature. It should be observed that the 
concepts of circulatory and noncirculatory flow components as developed 
in references 22 and 33 appear to have little meaning for wings with 
supersonic edges. Nevertheless, for convenience, these concepts have 
been utilized in the present method for wings with supersonic edges since 
it is believed that inclusion of the appropriate section lift- curve slopes 
and aerodynamic centers represents the principal aerodynamic effects on 
the calculated flutter speed of wings with supersonic edges. 

The section lift P and pitching moment Ma which are used in the 
present analysis may be obtained from similar expressions in reference 22 
by introducing variable section lift-curve slope Ci and variable 

aerodynamic center ac n . The procedure for making this generalization 
is as follows: 

First, the expressions for P and Ma used in reference 22 are 
written in the form 

l Noncirculatory 


P = -*pb c 


h + v n 9 + v n a tan - bale + v n T tan A^q. 


( 5 


1 


2jipv n bCQ 


Circulatory 


(3) 
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and 


Ma = + a2^e + v n f tan + *pv n b 2 (h + v n c tan A^j 

Kpb^a^h + v n 0 tan + npv n 2 b 2 ^e - abT tan - 


2«pv n b 


1 ' a + i)C 


> Noncirculatory 


Circulatory (i+) 


where Q is the downwash expression defined by equation (5a). These 
equations are, of course, based on the assumption that flow with small 
disturbances exists. 


Circulatory components . - Only the circulatory components of these 
expressions are changed. In the circulatory components of equations (5) 
and ( 4 ) the value 2 n for section lift-curve slope is replaced by the 
variable C7 * and the quarter-chord aerodynamic -center position 

^ ac n = ” 2^ re Pbaced by the variable ac n . The downwash expression Q 

must also be altered to include the effects of variable section lift-curve 
slope and aerodynamic center ac n . 

The treatments of the circulatory components of lift and pitching 
moment in references 22 and 33 are based on classical two-dimensional 
incompressible thin-airfoil theory, which indicates a section lift-curve 
slope of 2k and an aerodynamic center located at the quarter-chord posi- 
tion. The circulation strength is therefore related to the downwash veloc- 
ity at the three -quarter- chord position. This downwash as given in ref- 
erence 22 is 


Q = h + v n 8 + v n a tan Aea + b 


H( 


9 + v n r tan Aea 


(5a) 


and the distance between the bound vortex ( quarter-chord) the point 

at which the downwash boundary condition is applied (three-quarter-chord) 

C i 

For arbitrary Cj a ^ n , this distance becomes — b. (See ref. 2 k 


is b. 


2k 


for a detailed discussion of the application of the downwash boundary con- 
dition when bs other th an 2k . ) Then, if ac n (location of bound 
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vortex) is also arbitrary, the downwash condition is applied at the posi- 

( C l \ 

tion b( — + ac n J measured positive rearward from the midchord. (See 

fig. 1.) Then, in the expression for Q, the distance from the elastic 
axis to the point of application of the downwash condition b^i - aj is 


replaced in the present analysis by b 
present method. 


\ 2k 


+ ac n - a 


Then for the 


« 

Q = h + v n 9 + v n a tan Aea + b 


^ * a,n + ac n - a^6 + v n T tan Aea) (5b) 


Koncirculatory components . - The none ir dilatory flow components con- 
tribute to the lift and moment only a virtual mass effect which is com- 
paratively very small except at high frequencies. Since, as mentioned 
previously, the present method should probably be applied only to cases 
involving low to moderate reduced frequencies, it appears that the non- 
circulatory flow terms will constitute only a small fraction of the over- 
all section lift and moment. Now, the none ir dilatory components of sec- 
tion lift P and moment Mq, which are used in references 22 and 33 and 
shown in equations ( 3 ) and ( 4 ) of the present report are derived from 
the velocity potentials for unsteady two-dimensional incompressible flow 
about a flat plate. The virtual mass effects resulting from these non- 
cir dilatory flows are dependent only upon the velocity perpendicular to 
the wing surface and do not depend on the stream velocity as such. For 
low to moderate frequencies, the velocity perpendicular to the wing sur- 
face will be small compared to free-stream velocity. Therefore, for wings 
with all edges subsonic, any effects of compressibility on the magnitudes 
of the noncirculatory flow terms should be small, and the consequent 
effects on the section lift and moment should be of second order. It is 
concluded that, for wings with all edges subsonic, use of the noncircu- 
latory components of lift and moment in essentially the two-dimensional 
incompressible form should result in negligible error in the calculated 
flutter speed. 

In view of the relatively small magnitude of the noncirculatory flow 
components, the two-dimensional incompressible form is also used as a 
first approximation to virtual mass effects for wings with supersonic 
edges as well as for wings with all edges subsonic. At low reduced fre- 
quencies, the noncirculatory terms might even be completely neglected 
without introducing major errors into the calculated flutter results. 
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The section lift and. pitching moment used throughout the present 
investigation are made up of circulatory components generalized as pre- 
viously described and noncirculatory components used in the unaltered 
two-dimensional incompressible forms shown in equations ( 3 ) and (4). 

The resulting expressions are for the section lift 


P = -*pb 2 jji + v n 8 + v n cr tan Ag a - ba^9 + v n T tan A ea 
a,n pVnbCQ 

and for the pitching moment about the elastic axis 

Mq, = -«P b4 ^ + a 2 ^(9 + v n T tan Ag a J + *pb 2 v n (h + v n a tan Ag a j + 
wpb^afh + v n cr tan Ae a ) + itpb 2 v n 2 (e - abT tan Ae a ) - 



j Noncirculatory 
^Circulatory (6) 


/ Noncirculatory 


2jtpv n b 2 


I - (* - ””) c IT 


7 Circulatory (7) 


where the downwash expression Q is that defined in equation ( 5 b) . 

Note that in accordance with the discussion in reference 22 the terms of 
equations ( 3 ), (4), (6), and ( 7 ) associated with the variation of the 
velocity potential with lengthwise distance y ! are omitted. 

Substituting expressions (6) and ( 7 ) into the dynamical equations (l) 
and ( 2 ) and using equation ( 5 b) , together with the assumption of harmonic 
motion, yield two homogeneous flutter equations in the two unknowns h 
and 8_. The flutter determinant resulting from these flutter equations, 
expressions for the elements of the determinant, and the method used in 
solving the determinant for the flutter condition are given in appendix A. 
The remainder of the description of the method is concerned with the evalu- 
ation of the static aerodynamic parameters Cj a n and ac n and the cir- 
culation functions F and G which appear in the expressions for the 
determinant elements . 


Static Aerodynamic Parameters 

All calculations of static aerodynamic parameters are made by con- 
sidering the wing to be rigid and flat. 
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For M = 0 (reference).- The reference flutter speed is found 

for each wing by using n = and ac n = - ^ at M = 0. With 

these values the flutter equations (A12) and (A13) reduce to those given 
in reference 22. 

For 0 < M < 1.- At subsonic (and incompressible) speeds the span- 
wise distribution of C^ a is found by the lifting-line method of ref- 
erence 2b . In reference 2b charts of the necessary influence coeff icients, 
which facilitate rapid calculation of the loading, are presented. Although 
this method involves the application of boundary conditions and the eval- 
uation of load intensity at only seven spanwise stations, the resulting 
accuracy is considered adequate for present purposes, and the method is 
used because of its simplicity. Simple sweep theory is used to relate 


Thus, 


COS J ^ 03 < 


For all subsonic speeds the 


aerodynamic center is taken at the quarter-chord position 


( aCn = ■ i)- 


However, at subsonic speeds higher than those calculated herein it may 
become necessary to take aerodynamic -center changes into account. Details 
of the loading calculations are given in appendix B. 


For M > 1 .- At supersonic speeds when the wing leading edge is 
swept behind the leading-edge root Mach line (subsonic leading edge), the 
equations of reference 25 are used to calculate the static distributions 
of Ci a and ac n . The method of reference 25 is based on a superposi- 
tion of conical flows, and relatively simple formulas are given for cal- 
culating the loading. When the leading edge lies ahead of the leading- 
edge root Mach line (supersonic leading edge), the equations of 
reference 26 are used. Reference 26 is also based on , conical-flow 
concepts. These equations for lifting pressure have been used in 
integrals which yield section-lift and pitching -moment coefficients Cj 

and C ma ( and hence ac). The resulting expressions and details of 
their application are given in appendix B. The equations for and 

Cma given in appendix B make it unnecessary to refer to references 25 
and 26 for present purposes . 
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Circulation Functions 
The complex circulation function 


C = C(M,knr) = F(M,knr) + iG(M,knr) 


appearing in equations (6) and ( 7 ) and in the expressions of appendix A, 
modifies the otherwise-static circulatory components of lift and pitching 
moment to account for the effect of oscillation. The F function modi- 
fies the load component which is in phase with angle of attack, and the 
G function introduces out-of-phase load components. Values of the F 
and G functions used in reference 22 were those developed by Theodor sen 
(ref. 33) for two-dimensional incompressible flow about an oscillating 
airfoil. In the present investigation these values are again used for 
M = 0, but the functions must be modified to account for compressibility 
effects at M > 0. The modification used herein is based on loading 
functions for two-dimensional subsonic or supersonic flow about an oscil- 
lating airfoil as given by Jordan in reference 27- The relations between 
these loading functions and the F and G circulation functions are 
derived in appendix B. Although the flutter calculation is based on a 
consideration of sections normal to the elastic axis, the governing Mach 
number for the determination of the circulation functions is taken to be 
that normal to the leading edge. This choice of governing Mach number 
arises from the fact that the nature of the flow over a section of wing 
is influenced by whether the leading edge is subsonic or supersonic. 


Although it would seem straightforward to use the appropriate F q 
and Gq functions directly in the flutter calculations, this procedure 


gives poor results in comparison with experiment. (See figs. 3 and 9 > 


for example . ) The large phase angles 




of the complex circula- 


tion functions associated with two-dimensional compressible flow were 
found to be inappropriate for three-dimensional wings. It was antici- 
pated that if phase angles remained moderately small (i.e., if G remained 
fairly small relative to F)l, the calculated flutter speed would be rel- 
atively insensitive to changes in the magnitude of G. That is, if G 
is not large relative to F, the actual value of G is unimportant. The 


4he assumption of small phase angles implies an upper bound on the 
values of reduced frequency k^r for which the present method can be 
used. However, as previously mentioned the use of statically based load 
distributions also restricts the method to moderately small frequency 
values, so the present assumption imposes no further limitation. 
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predominant effect on the loading of changing Mach number would then lie 
in changing the magnitude of the in-phase component associated with F. 
The form of the complex function C which is used in the present calcu- 
lations is therefore taken to be 


C = 



F + iG = 



This function contains an in-phase component which is the same as that 
derived from reference 27 for two-dimensional compressible flow, but the 
associated phase angle is independent of Mach number. Hence, the phase 
angle is the same as that given by Theodorsen in reference 33* 

In order to investigate the validity of this reasoning some calcu- 
lations were also made by using 


C = 



= F c + iO 


Also, to investigate the sensitivity of the flutter calculations to dif- 
ferent forms of circulation- function representation, some calculations 
at the higher Mach numbers were made by using 


C = 




This function has zero phase angle, and its amplitude is the ratio of 
the magnitudes of the resultant vectors for compressible and incompres- 
sible flow. 

Further details of the circulation-function calculation are given 
in appendix B. The method for solving the final flutter determinant is 
given in appendix A. 
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RESULTS AND DISCUSSION 


Presentation of Results 

Flutter characteristics have been calculated by the present method 
(using three vibration modes) for wings with sweep angles from 0° to 
52 . 5°j aspect ratios from 2.4 to 7-4, taper ratios of 0.6 and 1.0, and 
center-of-gravity positions between 34 percent chord and 59 percent chord. 
The plan forms of these wings are shown in figure 2. The calculated 
results are compared with experimental data obtained in the Langley 
26 -inch transonic blowdown tunnel (refs. 28 to 31) and. in the Langley 
9- by 12-inch supersonic blowdown tunnel (ref. 32 )- 

Unless otherwise indicated the subsequent discussion deals entirely 
with calculated results obtained by using the complex circulation function 



Wing designation .- The three-digit system used to identify the wings 
with taper ratio of 0.6 is the same as that used in reference 30. The 
first digit in this system is the aspect ratio of the full wing to the 
nearest integer. The second and third digits give the quarter-chord 
sweep angle to the nearest degree. For example, wing 445 has an aspect 
ratio of 4, a sweep angle of 45°, and a full-wing taper ratio of 0.6. 
Since some of the wings discussed in this paper have identical plan forms 
but different center-of-gravity positions (ref. 31) > a single letter is 
appended to the plan-form designation to signify a shifted center of 
gravity. For example, wing 445 has a center of gravity at approximately 
46 percent chord, whereas the center of gravity of wing 445F is at about 
34 percent chord, and that of wing 445R is at about 58 percent chord. 

Wing 400 has a center of gravity at approximately 45 percent chord, but 
wing 400R has a center of gravity at about 59 percent chord. 

For the wings with taper ratio of 1.0, the same system is used, 
except that a fourth digit 1 is added to distinguish the taper ratio. 

For example, wing 4451 has a full-wing aspect ratio of 4, a sweep angle 
of 45°, and a taper ratio of 1.0. 

Flutter characteristics.- Calculated flutter characteristics v/Vr, 
aj/cqx, and k nr and the associated values of Vr, M, o^, and p are 

given in table I for several wings (see fig. 2) at several Mach numbers. 
The calculated values of V/Vr and co/(£& are compared with experimental 
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data in figures 3 to 14 and 15 to 2 6, respectively. The experimental 
flutter points shown were obtained at various values of density pj 
whereas, for a particular wing, all of the points calculated by the pres- 
ent method were obtained at a constant value of p which represented 
approximately an average of the experimental densities. For each experi- 
mental point, however, the normalizing Vp was calculated by using the 


appropriate experimental density. On the basis of previous experience, 
it is believed that normalizing the experimental flutter speeds in this 
manner essentially accounts for density effects so that the resulting 
(v/V R )exp is considered to be nearly independent of p, at least over 

the range of density variation which occurs herein. 


The static distributions of Ct „ and ac n used in obtaining the 

U/^ XI 

calculated flutter characteristics are shown in figures 27 to 35. For 
all of the flutter calculations presented in this report, the flutter 
modes of the wings were represented by a combination of the first torsion 
mode shape and first and second bending mode shapes of a uniform cantilever 
beam. 


The reference flutter speeds Vr used in references 28, 3 0 > 32 

for wings 430, 245, 400, 4001, and 7001 were calculated by employing 
only two degrees of freedom (first bending and first torsion). Since 
three-degree-of-freedom calculations yield values of Vr which are 

slightly lower than the two-degree-of -freedom values, the experimental 
V/Vr values for these wings have been multiplied by the ratio 

Vr for two degrees of freedom ^ ^ ^ caloiaated experimental 
V R for three degrees of freedom 

flutter-speed ratios as presented herein are normalized by Vr for three 
degrees of freedom. 


Flutter Speeds 

As shown in figures 3 to Ik, the flutter speeds calculated by the 
present method for all wings demonstrate a characteristic decrease as 
Mach number increases from 0 to near 1.0. This decrease is the result 
of increasing Cz a which is caused by compressibility at high subsonic 

speeds. It should be noted that at M = 0 the differences between the 
V/Vr values shown and the value 1.0 result solely from the effect of 

finite aspect ratio. As Mach number increases above 1.0, decreasing 
Ci and rearward shifting ac cause a rapid rise in the flutter speed. 

In the immediate vicinity of M = 1.0 the flutter-speed curves are shown 
dashed to indicate that this region is inaccessible to the present 
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calculations. This inaccessibility results from the breakdown of both 
subsonic and supersonic three-dimensional steady-flow wing theories near 
M = 1.0. It should be noted that the minimum value of V/Vr will gen- 
erally occur within this inaccessible region, and, hence, ( — ] can- 

^Yr Jm±n 

not usually be calculated by use of theoretical static aerodynamic coef- 
ficients obtained from the wing theories employed herein. It is possible 
however, to fair a reasonable curve through the neighborhood of M = 1.0 
by making use of the adjacent subsonic and supersonic calculated points. 
The extent shown for the dashed portion of the curves should not, of 
course, be interpreted as representing the limits of the inaccessible 
region. No attempt has been made to evaluate these limits, and the range 
shown in the figures is only illustrative. 

For all of the swept wings the calculated flutter-speed curves of 
figures 3 to l4 are in very good agreement with the experimental data at 
all Mach numbers. In general, the calculated curves actually lie within 
the scatter of the experimental data. For wing 445 (fig* 3) there are 
no experimental data in the range 1.4 < M < 1.75* However, the leveling 
off tendency demonstrated by the calculated flutter- speed curve in this 
Mach number range is in qualitative agreement with data for other similar 
wings . 

Comparison of the flutter-speed curves for wings 445> 445F, and 445R 
(figs. 3; 4, and 5) shows that the rather large differences between the 
center-of-gravity positions for these wings cause only very slight dif- 
ferences in V/Vr at subsonic speeds. At supersonic Mach numbers, how- 
ever, the data show that the characteristic rise of flutter speed with 
increasing Mach number becomes more rapid as the center of gravity is 
moved progressively forward. This behavior is also predicted by the 
calculated curves . 

The close agreement between calculated and experimental flutter 
speeds for wing 245 (fig* 6) is rather surprising in view of the small 
aspect ratio of this wing. In general, the use of a strip-theory type 
of analysis and uncoupled vibration modes for a wing of such small aspect 
ratio (panel aspect ratio = 0-91) would be open to question. The agree- 
ment in the present case may, therefore, be fortuitous. 

For most of the wings shown in this report no tip correction was 
applied to ac n to account for the forward shift of aerodynamic center 

within the tip Mach cone. (See discussion of tip corrections in appen- 
dix B.) For wing 430, however, the tip Mach cone covered so large a 
portion of the wing that it was considered necessary to apply a tip cor- 
rection to ac n . (See figs. 28(d) and (e).) At M = 1.15470, this cor- 
rection appears to be rather large. However, a preliminary calculation 
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at this Mach number without the correction to acn yielded a value of 
flutter speed only 13 percent higher than that shown. It appears, there- 
fore, that un less the tip Mach cone covers a large portion of the wing, 
the application of a tip correction to ac n is not necessary. 

For the lov-aspect-ratio unswept wings (figs, 10, 11, and 13) agree- 
ment between calculated and experimental flutter speeds is not as good 
as for the swept wings. For wing 400 (fig- 10) the agreement is fair up 
to about M = 1-0, but the calculated values overpredict the flutter 

speed by as much as 2— times at H -f- The magnitude of this error is 

believed to be related to the proximity of the local aerodynamic centers 
to the local centers of gravity and the fact that linear theory predicts 
an aerodynamic center that is too far rearward. This hypothesis is sup- 
ported by the results obtained for wing 400 with its center of gravity 
shifted from about 45 percent chord to about 59 percent chord (wing 400R) . 
Figure 11 shows that for wing 400R at supersonic speeds the calculated 
curve overpredicts the mean experimental values by only about 13 percent. 
The erroneous results obtained for wing 400 should probably not be inter- 
preted as indicating a limitation on the present method of flutter calcu- 
lation. Rather, these errors appear to arise from the well-known limita- 
tions on the use of linearized flow theory to calculate load distributions 
on wings of finite thickness. Wing 400 at supersonic speeds seems to con- 
stitute a very sensitive case in which a small inaccuracy in the location 
of the aerodynamic center leads to large errors in calculated flutter 
speed. In the case of wing 4001 (fig. 13) the calculated and experimental 
values are in very good agreement up to about M = 1.0. At supersonic 
speeds, where the local aerodynamic centers are shifted rearward toward 
the local centers of gravity, the theory again overpredicts the experi- 
mental values, this time by up to 37 percent. This deviation is not 
surprising in view of the fact that wing 4001 is not greatly different 
from wing 400. 

The calculated flutter speeds for the high-aspect-ratio unswept wing 
(wing 7001, fig. 14) are in good agreement with experiment throughout the 
Mach n um ber range . The improved agreement for this wing as compared with 
that for the low-aspect-ratio unswept wings may be caused to some extent 
by the decreased thickness of wing 7001 near the tip. Wing 7*^01 was 
tapered in thickness from 4 percent at the root to 2 percent at the tip, 
whereas wings 400, 400R, and 4001 were of constant 4-percent thickness. 

The flutter-speed curves shown in figures 3 "to 14 were calculated 
by using the complex circulation function 



20 


NACA m L57L10 




as described previously. The few points obtained by using the function 


C — F c + iO 


differ from the curves by no more than 7 percent. This close agreement 

supports the previously stated contention that if phase angles |tan~l £• 

1 \ F . 
are moderately small, the calculated flutter speed will be relatively 

insensitive to changes in G. Figures 3 to 14 also show that flutter 

speeds at the higher Mach numbers calculated by using the function 


\fc + G C f 
/=== ( 1 + 
\^I 2 + G! 2 



differ from the curves by no more than 10 or 11 percent. Although, as 
expected, the points calculated in this manner do not agree with experi- 
ment as well as the curves ^obtained with C = ps^Fj + iGj » y the small 

differences between them do point out the relative insensitivity of the 
calculated flutter speed to the form of circulation-function representa- 
tion used. 


In making the flutter calculations presented herein it was observed 
that for all but the highest subsonic speeds the circulation functions F c 

and Gq are not greatly different from the functions Fj and Gj of 

Theodorsen. At M = O .75 for the wings shown in figure 2, the use of 


C = F I + iGj 


instead of 


C = 



+ iG z ) 


changes the flutter speed by 


only about 4 percent or less. It would seem, therefore, that the modified 
circulation functions need be employed only at high subsonic and super- 
sonic speeds . 


■*-It should be clearly understood that the quantity tan - -*- — is the 

F 

phase angle of the complex circulation function C = F + iG and sho uld 
not be confused with any phase angles associated with the wing 
displacements . 
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Flutter Frequencies 


The calculated curves ^ith C = — ^Fj + iGj JJ of flutter frequency 

(figs. 15 to 26) indicate that for all of the swept wings the frequency 
is well predicted at subsonic speeds. At supersonic speeds the usual 
rise in frequency is predicted by the theory, but it occurs at Mach num- 
bers higher than those indicated by the test results. In general, the 
agreement between calculated and experimental, flutter frequencies is not 
as good as the agreement between calculated and experimental flutter speeds. 
The frequencies calculated for the swept wings by using 


' F c 2 + g c 2 ( 

If/ + G x 2 


1 + 10 ) 


are all excessively high, except at Mach numbers where the leading edge 
is supersonic or nearly so. At these higher Mach numbers the frequencies 
thus obtained are generally in better agreement with the experimental 


values than are the values obtained with 


v 


For unswept wing 4001 (fig. 25), the number of calculated points is 
not sufficient to indicate whether the pronounced dip in frequency, which 
occurs at high subsonic Mach numbers, is predicted by the theory. At low 

. „ , . . , . / ~ F C/„ . \\ 


supersonic speeds the calculated curves ^with C = —\F X + iGj JJ over- 

predict flutter frequencies by a substantial amount. However, the dif- 
ferences between theory and experiment become much smaller at the higher 
supersonic speeds, except in the case of wing 400 (fig. 22). The fre- 
quencies as well as the flutter speeds of wing 400 are overpredicted by 

a factor of nearly 2—. As in the case of the swept wings the frequencies 

2 

for the unswept wings obtained by using 


I p 2 . n 2 

f *c + °c , 


;i + io) 


Pi 2 + G i 2 


are all excessively high. 
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Limitations of the Method 

Although the limitations of this method have not been fully evalu- 
ated, some of the more important restrictions may be qualitatively 
discussed. 

Frequency range . - As stated previously, the use of spanvise load 
distributions based on lift-curve slopes and aerodynamic centers calcu- 
lated from steady-flow wing theory imposes an upper bound on reduced- 
frequency values for which the method can reasonably be used. No attempt 
has been made to determine the upper limits of reduced frequency for which 
the method is usable, but good results for values of knr up to 0.2 are 
shown herein. 

Mach number range . - The nature of the equations for the circulation 
functions (eqs. ( B 3 8 ) and (B39) or (BkO) and (B4l) ) shows that at M j.k = 1, 
the circulation functions become Fq = Gq = 0. This implies that a small 
range of Mach number in the immediate vicinity of M-p.^ =1 is inaccessible 

to the present method. This is not a serious limitation, however, because 
a curve of flutter speed or frequency can be reasonably faired through 
this inaccessible region by making use of adjacent points. For the wings 
calculated in this report, there appear to be no sudden or extreme fluc- 
tuations of flutter speed or frequency in this region. 

The limitations on Mach number range appertaining to the particular 
steady-flow wing theories used are, of course, carried over to the flutter 
calculation. In general, this carried-over restriction will exclude free- 
stream Mach numbers in the immediate vicinity of 1.0, as was mentioned 
previously. 

Flutter modes .- The use of uncoupled modes in combination with a 
strip theory involving strips normal to the elastic axis is not an 
essential requirement of the present method of flutter calculation. An 
analogous calculation procedure would result from the use of coupled 
modes together with streamwise strips. Flutter modes which involve 
significant amounts of camber deformation obviously cannot be treated by 
the method in its present form. As mentioned previously, all flutter 
calculations presented herein were made by using the mode shapes of a 
uniform cantilever beam. Since the results of the flutter analysis are 
not very sensitive to slight changes in mode shape, such a procedure is 
reasonable as long as aspect ratio and especially taper ratio are not 
too small. 

Plan-form range .- The strip-theory concepts which are employed in 
the present method also impose plan-form limitations. When aspect ratio 
or taper ratio or both become so small that the variables (notably herein, 
aerodynamic loading and circulation function^ «-sociated with a given 



NACA RM L57L10 


25 


section of the wing cannot be treated by strip theory, then the present 
method is no longer usable. 

Center- of -gravity position .- Although the influence of different 
center-of-gravity positions was investigated for only two plan forms 
(wings 4^5 and 400), it appears that, in cases for which the local centers 
of gravity are located close to the local aerodynamic centers, linearized 
flow theory should be employed only with great caution. This limitation 
is not peculiar to the present method. It would apply to any flutter cal- 
culation for which the aerodynamic loadings are obtained from linear 
theory. 

At subsonic speeds, neither the swept nor the unswept wings demon- 
strate any appreciable sensitivity of V/Vr to center-of-gravity posi- 
tion. This result would be expected since at subsonic speeds local aero- 
dynamic centers are at or near the quarter- chord position and are not in 
proximity to the local centers of gravity. 


CONCLUDING REMARKS 


A method has been developed for calculating flutter characteristics 
of finite-span swept or unswept wings at subsonic and supersonic speeds. 
The method is basically a Rayleigh type analysis and is illustrated herein 
with uncoupled vibration modes although coupled modes can be used. The 
aerodynamic loadings are based on distributions of section lift-curve 
slope and local aerodynamic centers calculated from three-dimensional 
steady-flow theory. These distributions are used in conjunction with 
the "effective 1 ' angle -of -attack distribution resulting from each of the 
assumed vibration modes in order to obtain values of section lift and 
pitching moment. Circulation functions modified on the basis of loadings 
for two-dimensional airfoils oscillating in a compressible flow are 
employed to account for the effects of oscillatory motion on the magni- 
tudes and phase angles of the lift and moment vectors. 

Calculation of subsonic and supersonic flutter characteristics for 
12 wings of varying sweep angle, aspect ratio, taper ratio, and center- 
of-gravity position and comparison of the results with experimental 
flutter data indicate that the present method gives generally good flutter 
results for a wide variety of wings. The method is, however, subject to 
the following limitations: 

CD it is probably not applicable at high values of reduced frequency, 
although good results are shown for values of reduced frequency up to 
about 0.2. 
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(2) It cannot be used at free-stream Mach numbers in the immediate 
vicinity of 1.0 nor in the immediate vicinity where the Mach number com- 
ponent normal to the leading edge is 1.0. However, flutter speeds and 
frequencies may be interpolated through these regions. 

(3) The use of a strip-theory approach and the absence of camber 
flexibility preclude treatment of wings with low aspect ratio and low 
taper ratio (e.g., delta wings). Good results have been obtained, how- 
ever, for a 1+5° swept wing with a panel aspect ratio of 0-91. 

(4) Caution must be used when applying the method to wings for which 
the local aerodynamic centers are close to the local centers of gravity . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., November 26, 1957* 
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APPENDIX A 


DERIVATION OF THE FLUTTER EQUATIONS AND FLUTTER DETERMINANT 


Flutter Equations 

Basic assumptions .- The dynamical equations used in the present 
method are essentially the same as those derived in reference 22, except 
for changes in the expressions for lift P, pitching -moment and 

circulation functions F and G. The general assumptions appertaining 
to the method of reference 22 thus apply herein also. Briefly, the 
assumptions made with regard to the equations of motion are as follows : 

(1) The elastic axis of the wing is approximately straight and 
the oscillatory motion may be represented by a combination of the 
uncoupled bending and twisting vibration modes of the wing with respect 
to this elastic axis. 

(2) The wing root is treated as though it were clamped along a 
line normal to the elastic axis and passing through the intersection 
of the elastic axis and the root chord. 

(3) The analysis is based on geometric, structural, and aerodynamic 
quantities associated with sections normal to the elastic axis. These 
assumptions are discussed in detail in reference 22. 

Application of Lagrange's equations .- The dynamical equations 
result from the application of Lagrange ' s equations of motion to the 
flutter problem. For simplicity, the flutter equations are derived 
herein for the case of one bending mode and one torsion mode. 
Generalization to an arbitrary number of modes is easily accomplished 
in the flutter determinant as will be illustrated. (The notation of 
of ref. 22 has been followed where possible.) In the present method 
the appropriate expressions for kinetic energy 


T = 



1 2 


dy' 



2 dy' + 



(Vy’)j d y' 


(Al) 
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potential energy 


u * iV* 




raf h 2 (iy * 


+ 2 0 2 
2 cl — 


r 
J 0 


Ve 2 ^ ' 


and virtual work 


5W = Q h &h + Q 0 60 

are the same as those of reference 22. The generalized forces are 
in the form 


% ■ f o \r - W f hi) f h dy ' 

and 

«e - f ’ (t - W I 6 

Substituting these expressions into Lagrange's equations 


d /8 t\ dT au _ 
dt^hj - a» ah - ^ 

and 


d /^t\ dT dU _ Q 

dt^e j (30 c30 e 

and assuming harmonic oscillations 


h = 


f h (y')]h = [fh(y')]v laJt 


(A2) 


(A3) 

left 


(A4) 


(A5) 


(A6) 


(at) 


(A8) 
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and 


0 = [ f 0 (y ')] 0 ■ [ f e(y'j] e c 

lead to the equations of flutter 


a iojt 


(A9) 


-f -( 1 +iig h ) - i 


L 




-b. 


r 1 i/jlV 

*' 0 K ( b rj 


Vh f e d y 


1 • ^ /o' • 0 


(A10) 


and 



i 

*pb r 3a£ 


£ W - 0 


(All) 


In the calculations of the present report, uncoupled beam bending and 
torsional mode shapes hj. and a.j are used for the flutter deflection 

functions f^ and f g . The introduction of uncoupled modes into the 
flutter equations is discussed in detail in references 22 and 

Expressions for the elements of the flutter determinant resulting 
from equations of the type (A10) and (All) are given in the following 
section both for the case of an arbitrary number of vibration modes 
and for the case of one torsion and two bending vibration modes as used 
in the present analysis . 


The Flutter Determinant 


Inserting equations (6) and (7) into equations (A10) and (All) 
and using equations (A8), (A9) , and (5b) yield two homogeneous equa- 
tions in the two unknowns h and £, which may be written in the 

form 


Ah + B9 = 0 
Dh + E 9 = 0 


(A12) 


and for a nontrivial solution to exist, 


A B 
D E 


0 


(A13) 
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Now if u bending modes and v torsion modes are employed, the 
elements of the flutter determinant (A13) will become matrices k^y 

Bij, D^j , and E i j , such that 


ij 

B U 


E U 


= 0 


(A14) 


The solution of equation (Al4) gives the conditions of flutter 
(flutter speed and frequency) . The procedure for solving this determi- 
nant is given at the end of this appendix. Expressions for typical 
elements in the matrices A^j, B^j, ^ij* and ^ij are as f°H- ows • 


l ii 


(^ 1) 2 (1 + ig h \ " 1 l f 7 B^ 4 . 2 dn - 1 f B 2 h ± 2 &T\ 
\ “ ) \ h i/ \ Jo K J 0 


, _ r 1 o b r tan Ao a n „ dh i . , 

1 - rr— / c i Bh, 2 dq + 5 C l B h i dT > 

* Vr J 0 l a,n 1 * 0 a ’ n dn 


+ 


b r tan 


Aen P d p 

-2* / B d -r-i h.dq 
J n dq i 


^nr J 0 

(i = 1, 2, 3, . . . u) 


knr J 0 


AiJ --if 1 *Vj*< * 1 1 £ f Q C l a ,„ Eh i h J 4 '' 


b„ tan A~ a p P 1 _ dhj 

_£ — / C 7 B - — i h.dq + i 

* k^ 2 Jo a > n dT l 1 


( 


>j* tan Aea P 

‘ B 2 ^ 

^nr 0 

dq 

i = 1, 2, 3, 

. . . u 

j = 1, 2, 3, 

. . . u 

i / J 



h^ 
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B iJ = 



^ B^ajdri 


+ v_c_ r 1 

* 0 


C, Bh^cudTj + 


v ± r* 

* Vr J o 


B* 


^a,n 


2* 


+ ac n - a h. 




v 


j B^ahja.drj 

Jo J 


b r Z 


Kr J 0 


dT] 


b r tan A^g. q 

* w 


- f Cj E 
! J o a-,n 


+ ac . 

2 jc 


Wj 

- - a jd^T h i dT > - 


i 


b r "tan A^a 
k nr 





da. 


/ i = 1, 2, 3, . . . u\ 

U = 1, 2, 3, • • • y) 



b r 2 tan Ag a 

Jt 




dh 1 

a M i^r a±dT1 " 


i 


b r tan Ae & 
^nr 



B 5 a 


dhj 

dTi 


a^dt) 


/i - 1, 2, 3, • • • v\ 
= 1, 2, 3, . . . uj 



50 


NA.CA RM L57L10 



(equation continued on page 31) 



NACA RM L5TL10 


31 


, 2 1 pl ,°l 


'■nr J o 


+ ac - a/a.a,dT| 


f-Jo + ac - - a )( a ■ a °^ ^ + 


b r ^tan 


> J 0 \ 2* n / I! i 




i * i, 2 , 3 , ... ^ 

j = 1, 2, 3, . . . v 

i 4 J ; 


In the special case of three degrees of freedom (first and second 
bending and first torsion modes) used throughout the present investiga- 
tion, the flutter determinant (eq. (Al4)) becomes 


•11 

A 12 

B 11 



‘21 

A 22 

B 21 

= 0 

(A15) 

>11 

D 12 

E 11 




The elements of this determinant can be conveniently expressed in the 
forms 


g , „ r\ x _L 


An - a, + a 2 — + a 5 — + i - 


+ a 5 — - Rtf 

K nr / 


A-.o = fa + b 5 + b, -^o) + if}— " b 2 ^ + b 3 kr G 


^12 = ri * U 2 ^ - 3 knr 2 
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J 12 
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G 
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m. 





where 

C = F + iG 


and 


b r o) 

Kr = v^~ = v cos Ae a 
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„ _ b r tan 

a 5 * 


© 


a 4 = b r tan Ag a © 


>1 = -i 




b r tan Ag a 

b* = 

3 if 


<u 


b]j_ = b r tan 


\a © 


=1 “ wk © + V 


-b„Z 


c„ = 


2 * 


c = hi (u) + b r tan ^a 

3 « ^ * 


c = b r l @ - b r 2 tan A^ 
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v 
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b„ tan 
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m c = b ''tan A 
5 r ea 




_±_Nf 

«pb r ^ \ °V 




and the circled numbers represent the following Integrals : 
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©- 

r 1 dh i 
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These integrals are easily evaluated numerically. Only about one-half 

of these integrals contain C ? or ac^ . Hence, only these integrals 

L a,n 11 

change with Mach number. For a given wing the remaining integrals may 
be evaluted once for all. Note that the integrals are independent of 
density p. The density appears only as a multiplying factor in a^, 

c x > e 1 > f ± > g-^ h 1> m i , R g , and R^. 


Solution of the Determinant 

For a given wing at a given Mach number the three-by-three flutter 
determinant (eq. (A15)) was solved for Z on an electronic digital 
computer for various values of the parameter k nr (and associated values 

of F and G) . This evaluation of Z yielded values of — — and g 
corresponding to the various values. A plot of g against 
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v 

then gave the value of k_ and — — for which g = 0. These values 

b r“k 

b-pU) 

*nr " — 
v n 

and 

v n _ Vn 
b r^ b r“a 

define a flutter point. Then the flutter speed V is 


V = 


' n 


b„co„ 
r a 


b i"\x 


cos 




and the flutter frequency a) is 


o> = ^nr x 


'n 


r a 


x 
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APPENDIX B 

PROCEDURE FOR MAKING FLUTTER CALCULATIONS 


The following procedure was used in making the calculations presented 
herein. 


Summary of Required Information 


First, a summary sheet is set up similar to that shown in table II. 
The entries on this sheet represent all the information necessary for the 
evaluation of integrals (T) to (27) > coefficients aq to mjp, and Rq, 
R 2 , and R3 listed in appendix A. These coefficients together with the 

circulation functions F and G (calculation of which is discussed at 
the end of this appendix) permit evaluation of the determinant ele- 
ments Aqq to Eqq and, hence, solution of the flutter determinant as 

described in appendix A. 


Columns (l) to (5) of the summary sheet contain wing mass and elastic 
parameters which, in the present case, were determined experimentally. 

All of the experimental flutter data shown herein were obtained with the 
wings mounted on a fuselage. (See refs. 30 to 3^--) The calculations 
were therefore made considering the wings to be cantilevered from the 


side of the fuselage which was assumed fixed. The quantities <%,, 




— — ) , and l listed at the top of the summary sheet are also measured 

values. Column (6) contains values of B = — , the nondimens ionalized 

b r 

semichord measured perpendicular to the elastic axis. The nondimens ion - 
alizing value b r is the semichord b at station tj = 0.75- Values 

of b may be obtained from the following equations: 


b = 6 TE ~ *LE (el) 

S ( K 1 + K 2 ) + ( K 1 " K 2 ) a 


< 
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where 


£le = tan Ale 


(B2) 


4 

|rjrg = — + T tan Apg 

A p (l + A p ) 


(B3) 


s = 1 cos Aea 


% = 


cos (Ale ~ Aea) 


cos Ai 


UE 


K 2 = 


C0S ( A ea - Ate) 


COS Ajrg 


(B4) 


tan Ay -ft = 


tan A c f^ + 


1 1 - A 
A 1 + A 


tan Aijrg 


tan Aj^e 


4 1 - A 
A 1 + A 


(B5) 


The geometrical quantities appearing in these equations sire shown in fig- 
ure 3 6 . Note that in equation (B3) the values of aspect ratio Ap and 

taper ratio Ap to be used are those obtained by considering the side of 

the fuselage to be a reflection plane. In equations (B5) it is immaterial 
whether A and A are obtained by considering the reflection plane to be 
at the side of the fuselage or at the fuselage center line. 

Columns ( 7 ) to (12) of the summary sheet (table II) are the ampli- 
tudes and slopes of the uncoupled vibration mode shapes. These mode 
shapes may be calculated for the particular wing by any of the methods 
given in references 36 and 37- However, since flutter speed is not highly 
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sensitive to slight changes in mode shape, the mode shapes for a uniform 
cantilever beam may be used if aspect ratio and taper ratio are not too 
small. All calculations in the present report were made by using the 
first torsion and first and second bending mode shapes for a uniform canti- 
lever beam as given in table III and figure 37- (Equations governing har- 
monic bending or torsional oscillations of a uniform cantilever beam are 
derived in ref. 37*) Table III contains all combinations of these mode 

shapes which are required for the calculation of integrals (T) to (2j) . 

Also presented in table III are the integrals of these mode-shape com- 
binations which are useful in evaluating the integrals for untapered 
wings . 

Columns (13) and (l4) of table II represent the distributions of 
static aerodynamic parameters at a given Mach number. 


Calculation of Static Aerodynamic Parameters C^ a n and ac n 

The values of local lift-curve slope C7 are obtained for sec- 

tions normal to the elastic axis by applying simple sweep theory to C^ 

<x 

values for streamwise sections. Thus, 




a > n cos Aea 


(B6) 


The use of simple sweep theory together with values of C for stream- 

wise sections results in C7_ _ values different from those obtained by 

direct integration of pressures over sections normal to the elastic axis. 
However, the resulting discrepancies are negligibly small except near the 
wing root where deflection amplitudes are small. (See fig. 27 (e).) The 
use of simple sweep theory should thus cause negligible errors in the 

values of the integrals © to (2j) . The local aerodynamic centers ac n 

in units of semichord b and measured perpendicular to the elastic axis 
are found from 
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ac n = 


ac 


^TE ~ ^LE 1 


b/ £ 


(l + a)K x 


cos Ae a + a 


ac 


(ki + K 2 ) + (ki - K 2 )aJ - (1 + a)Ki 


cos + a (B7) 


Distributions of Cj a n and. ac n for all the wings calculated are shown 
in figures 27 to 35* As indicated in figures 27(e) and (f) the values of 
C^ a and ac n used in the flutter calculations do not always lie on the 

curves of Cj n and ac n distribution. The integrals (l) to (2j) 

(appendix A) are evaluated numerically by using values of mass, elastic, 
and aerodynamic parameters at T) = 0.05 "to 0*95> in increments of 0.10. 

The required values of Cj a n and ac n therefore are average values 

over the T]-intervals 0 to 0.10, 0.10 to 0.20, . . . 0-90 to 1.00. These 
values do not coincide with the Cl a ^ n and ac n distribution curves 

near points of sharp change . 

Subsonic free stream .- In the case of subsonic free-stream velocity, 
the spanwise distribution of Cj a is found by the method of reference 24. 

For these subsonic loading calculations, the full wing is considered. 

That is, the reflection plane is considered to be at the fuselage center 
line, and the presence of the fuselage is neglected. The effect of the 
fuselage on the actual loading is felt primarily near the wing root. 

Since deflection amplitudes are small near the root, the overall effect 

of the fuselage on the integrals © to (§7) should be negligible. Since 

the loading distribution is computed for the full wing including fuselage 
intercept and since the distribution only over the wing panel is required 
in the flutter calculation, the full-wing distribution of Cj a is plotted, 

and values are read off at stations corresponding to q = 0.05, 0.15, 

. . . O.95 of the wing panel. (See fig. 38.) For subsonic free-stream 
velocity, ac n = -O.5 is used throughout. This value corresponds to the 
aerodynamic center at the quarter-chord of a section normal to the elastic 
axis . 


Supersonic free stream .- For supersonic free stream, the cases of 
subsonic leading edge and supersonic leading edge are considered. 
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(l) Subsonic leading edge: In the case of supersonic free-streain 

velocity and subsonic leading edge, the spanwise distributions of Cj a 

and ac are found by the method of reference 25. For these calculations, 
the wing is treated throughout as though the side of the fuselage is a 
reflection plane. This assumption seems reasonable since in the linearized 
theory of reference 25 the distribution of loading on the wing panel is 
dominantly affected by Mach waves emanating from the wing-fuselage 
juncture . 

When the leading edge is subsonic and the trailing edge is super- 
sonic, as in sketch 1, 



the expressions for streamwise C^ a and ac take a very simple form 


C 





4 

E tan Ajjj; 


pTE + ^LE 

V^TE ~ £LE 


(B8) 
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second kind. Expressions for ^ TF are 6^ ven by 


equations (B2), (B3), and (B5) . The numerical subscripts throughout 
refer to the loading areas in the appropriate sketch. Note that for 
this condition ac is a function only of wing geometry and that Mach 
number affects C^ a only through the function E. Equations (B8) 

and (B9) contain no provision for accounting for the loss of loading 
within the Mach cone from the tip leading edge. The procedure for 
applying tip corrections is discussed subsequently. 

When the leading edge and trailing edge are both subsonic, as in 
sketch 2, 



Sketch 2 
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the loadings indicated by equations (B8) and (B9) must be corrected to 
account for the loss of loading behind the root trailing-edge Mach line. 
For this condition 


C 7 = Ci + dC 7 , 

L a 4 a,l t a,3 


“ Cm a,l + > 


ac = — — 


-a 


(BIO) 


where and 


are obtained from equations (B8) and (B9), and 



where K = Kl \ /I - 


o 

tan ApE 


Is the complete elliptic integral of the 


first kind, Flcp, \ 1 - 


the first kind, and 


p 

tan Ape 


is the incomplete elliptic integral of 


4m 


+ pT) 

Ap(l + Ap) 


(B12) 



46 


NACA RM L57L10 


cp = sin"l 




Ap(l + Ap) 


(B13) 


Equations (Bll) represent only the "symmetric" trailing-edge correction 
discussed in reference 25- However, this quantity is considered suf- 
ficiently accurate for present purposes. The integrals in equations (Bll) 
are evaluated numerically. 


For T} stations near the wing tip the loadings given by equa- 
tions (B8), (B9) , or (BIO) must be corrected to account for the loss of 
loading within the Mach cone from the tip leading edge. When the leading 
edge is subsonic, as in sketch 3, 



Sketch 3 
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these corrections are made as follows: First, the spanwise locations of 

points P, S, and L (see sketch 3) are found from the equations 



Tls = n P + (^TE - Sea)n=Tj p sijl cos *** 


(B15) 


where 



(^ea - £le)tv= 1 “ (ijr^ 1 + a *l=l) K i 

(See fig. 36 .) The more inboard (measured parallel to the elastic axis) 
of the points S and L represents the tj station at which the tip 
effect first begins to be felt. 

The load intensity on the wing rises from trailing edge to leading 
edge and approaches infinity at the leading edge. Therefore, if < t^, 

the loss of loading caused by the tip will begin at the trailing edge 
where load intensity is relatively low and gradually extend forward into 
a region of high load intensity as the tip is approached. The loss of 
loading outboard of r\ Q will thus produce a curve of q as a func- 

tion of t) which has negative curvature as well as negative slope. (See 
fig. 39(a).) Now the static aerodynamic loading parameters are intro- 
duced into the flutter equations through strip theory which implies that 
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the loading has a quasi-two-dJmensional character. For a swept wing an 
exact stripwise evaluation of loading near the tip would therefore have 
questionable significance since neither the wing plan form nor the pres- 
sure distribution is quasi-two-dimensional in that region. In view also 
of the difficulty in performing an exact stripwise integration of loading 
near the tip, a reasonable fairing of the Cj a n curve is considered 

adequate, even though this fairing occurs at spanwise locations where 
wing deflection is greatest. For the case of qg < q^ (fig* 59(a)), 

the approximate curve used is geometrically derived from that obtained 
by streamwise integration of loading in the tip region. The geometrical 
derivation consists of applying a constant stretching factor to the curve 
obtained by streamwise integration of the tip loading in order to fit 
this curve to the known loading at q = qg. The appropriate equations 

for this streamwise calculation are equations ( 6 ), (l 5 )j and ( 26 b) of 
reference 25 . No reflections of Mach lines from plan- form edges are con- 
sidered. For wing 445, flutter speed determined by using this type of 
fairing and that obtained by using exact stripwise integration of tip 
loading differed by only 0.6 percent. 

If rjg = the curve of n has a sharp discontinuity at 

q = Tjg = q^. (See fig. 39(h).) In this case an accurate representation 

of the loading in the tip region can be obtained with the aid of figure 7 
of reference 25* This figure gives the loss of lift across the tip Mach 
line. A straight line is used, as in figure 39(h), to fair the Cj a n 

curve to zero at rjrp - The value of q ij< is given by 

n T = 1 + 6te “ ^ea) T]=l sin Aea cos ^a ( B1 7) 


If T]g > tj^, the region of high load intensity near the leading edge 
is lost first, so that the curve of C} a Q against r\ has a steep nega- 
tive slope just outboard of but has also a positive curvature (as in 

fig. 39(c)). In this case a straight line is used to fair the C n curve 
between tj^ and T) = 1. In no case is any loading outboard of tj = 1 
used in the flutter calculation. 

In general, no tip correction was applied to ac n since such cor- 
rections would occur in only a small region. For wing 430, however, the 
point was so far inboard that it was considered necessary to apply 

a tip correction to ac n . (See figs. 28(d) and (e).) This correction 



NACA RM L57L10 


^9 

was obtained in the same manner as the correction for Ct . That is, 

ex, n 

the correction was determined from streamwise integration of lift and 
pitching moment from which ac and hence ac n were found. 

(2) Supersonic leading edge: when the leading edge is supersonic, 

as in sketch 4, 



the spanwise distributions of Cj a and ac are found by the method of 

reference 2 6 . Again the wing is treated as though the side of the fuse- 
lage is a reflection plane. Values of and ac n are found from 

OC/, n 

Ci a and ac by applying simple sweep theory as described previously. 
The procedure for finding C^ a suod ac as First, find the 

spanwise locations of points P, 0, and Q (see sketch 4) by using the 
equations: 
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tan Aj-p + p - 


Tip - 


At.( 1 + Ap) 


V[ ± * A P 


tan Aipp + p 


no = 


tan A t ^ + p 

~ip 


no = 


^ Ap^l + Ap) p - tan Arpg 
If point 0 lies on the wing, then for 0 ^ rj ^ Tjp, 


C-, = 


l a | | 

\|P 2 - tan^Ay^ (|tE " ^Le) 


C m i ■ 5le ) + /, 


^TE /C p ,i 


6 Ml 


m 


and 


c % “ 


\f & 2 ~ tan 2 ^LE^TE ’ ^Le ) 2 
For Tip < T) < r\ 0 , 




$LE 


^TE " ^LE 


C *a 


C 7 = 


\j$ 2 - W^(l» - ^le) _ 


(IM, - Sle) * / ! " 2 + (*“ 

V V^ay J|m 2 \Cp, 


2D 


l<u 


and 


: 


\| P 2 - tan 2 ^^ - 




(B18) 


(B19) 


(B20) 


(B21) 


(B22) 


(B2J) 


(B24) 
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For no ^ ^ ^ IQ* 


C l<x = 


and 


\/p 2 - tan 2 ALE (^TE ' *Le) [ 


Ctei 


dll 


(B25) 


Cow ~ 


\j $ 2 ~ tan 2 AL£ (^TE ' ^Le ) 2 2 




2Dj 


<U 


S LE 


^TE " ^LE a 


C, (B26) 


For t]q = t] = 1> 


C *a = 


\jp 2 - tan^LE^TE “ ^Le) [ 


( tM2 ■ ii£ ) + c fei) 45 


(B27) 


and 


C ma 


\Jp 2 - tan 2 A LE (^TE " ^Le)* 

where, as before, 




Ile 


Ite - Ile a 


Cj (B28) 


^LE = n tan a LE 


> 


Ite = 


’f 1 + a p) 


+ r\ tan Arjrg 


(B29) 
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and 


Also, 


where 


Imi - Pn 

Im 2 = tan a le + P(i - tj) 


(B50) 


' Cp, 1 \ 1 i | 

= — cos _1 [R + 


St)‘ 


V C P,2D/ 11 


2 2 2 

| - taxi At^t) 


C P^2 ■; _ 1 / I - T \ 


, = — COS 1 

Cp,2Dy * U - tan 


(B31) 


'Cp,}\ C P) 1 C P; 2 

Cp,2Di "lCp,2Di \C V ,2D " 


R = — tan^Ajj; - 1 
p 2 


S = (R - l)tan 2 A] J E 


(B32) 


T = tan A Ti ^ + ( 2(3 + tan 


^2p + tan Aj^^l - i}) 


If point 0 lies behind the trailing edge, then rjp > and 
and Cn^ are obtained as follows: 

For 0 < rj ^ t]q, and equations (B21) and (B22); 

for t|q < n s n P , 


t 
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c l 


a 



k 

taia^A-LE 


(B33) 


Cma = - ■ ■ k x | (B34) 

UP 2 - tan^LE 

For Tjp 3 T) — 1 , and. Cn^ are given by equations (B27) and. (B 28 ). 

All of the integrals in the foregoing expressions for Ci a and 
are evaluated numerically. 

It should be noted that for the case of supersonic leading edge if 
T)S < no separate tip correction is necessary. Approximately correct 

values of n and ac n in the tip region (t] > Tjg) are obtained by 

applying simple sweep theory to the values of and ac resulting 

from equations (B23) to (B 28 ) . Loadings of this type are shown in 
figures 28(e) and 34(d). If T]g > t]-^ then the curve is faired 

with a straight line between and r\ = 1 . (See fig. 39(c).) In this 

latter case, equations (B23) to (B28) need not be evaluated. 


Circulation Functions 

As mentioned in the body of this paper, the circulation functions F 
and G, which appear in the determinant elements listed in appendix A, are 
obtained from aerodynamic coefficients given in reference 27 for two- 
dimensional airfoils oscillating in compressible flow. (Similar coeffi- 
cients for supersonic speeds only are also given in ref. 38*) These coef- 
ficients l a , l z , rrvx, m z are defined in reference 27 so that 


P = -2bp^(h m l z + 0Z a ) (B35) 


(B36) 


and 


Ma = (2b) 2 pv 2 ^h LE m z + 91%) 
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in the notation of the present paper, where hj.Tr. is the value of trans- 
lation deflection at the leading edge. Now another expression for lift 
in the case of two-dimensional compressible flow may be obtained from 
equation (6) by deleting the terms containing a and t. Thus, 


P = -«pb^(v9 + ii - ba9) 



+ ac n 



(B37) 


where 


= 

*a,n R 


> for M < 1 


ac n = - g 


and 


C *a,n p 

> for M > 1 

ac n = 0 

> 

Expressions for the circulation functions in terras of the aerodynamic 
coefficients of reference 29 may be obtained by equating expressions (B35) 
and (B37) . Equating the two expressions for P (eqs. (B35) and. (B37)), 
using 


9 = ioj9 
h - ioii 

0 = -o£e 
ii = -o£h 
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for simple harmonic oscillation, and noting that 



where 


i<x - U(^Anr) - 

l 2 = ^z(^>^nr) - ^z* + 

Considering only the pitching oscillation, that is, putting h = 0, 
permits simplification to 
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or 



Analogous expressions for F c and G c could be obtained by equat ing 
expressions for pitching moment instead of lift P. It was indi- 


cated previously in this report that use of the present method for pre- 
dicting flutter characteristics should probably be restricted to cases 
for which k^ is moderately small. Therefore, the knj.5 term in equa- 


tion (B39) may be dropped. 


Furthermore , 


the factor 



ac n 


does not vary greatly with Mach number except in the immediate vicinity 
of M = 1, and this vicinity is inaccessible to the present method. 
Therefore, since this factor is always multiplied by k^ or k x\x^~ > 

only small error will be introduced into the circulation functions by 


talking throughout 


'^a,n 


+ ac n - a I = — , which is the incompressible 


2 * 

* 

flow value with a = 0. The value a = 0 implies torsional oscillation 
about the midchord. Equations (B 38 ) and (B39) then reduce to 



(Bi4-0) 


CONFIDENTIAL 
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These expressions for Fq and. Gq are independent of wing parameters 

and depend only on Mach number M and reduced frequency knr* As men- 
tioned in the body of this report, when the two-dimensional circulation 
functions Fq and Gq are used in flutter calculations for three- 

dimensional wings, the functions are defined by the Mach number normal 
to the leading edge. Thus Cq becomes 

Cq = Cc^LEjknp^ = FQ^Mj^jk^j + iGc^M^jk^ 

A typical comparison of Fq and Gq calculated from equations (B40) 

and (B4l) with those obtained from equations (B 38 ) and (B39) is shown in 
figure 40 . Values of Fq and Gq were obtained from equations (B 38 ) 

and (B39) for two positions of aerodynamic center: ac n = 0 (the two- 

dimensional supersonic value) and acn = -O .325261 (the value at the 

station q = O .75 of wing 445 at M = 1-75) • The results in both cases 
closely approximate the results from equations (B40) and (B4l). The dif- 
ferences between the three sets of Fq and Gq curves shown in figure 40 

would result in less than 1 percent difference in the calculated flutter 
speed for wing 445 . Since calculated flutter speed is only moderately 
sensitive to small changes in the circulation function values (see fig. 3)> 
the circulation functions used throughout this investigation were cal- 
culated from the simplified equations (B40) and (B4l) . 

Some typical curves of Fq and Gq are shown in figure 4l, and 
the combinations ( 21 ^' - * z ') and ( 2 Z a " - l z ") used in equations (B40) 
and (B4l) are plotted in figures 42 and 43, respectively. 
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L Values listed are for uniform cantilever 
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(a) Relations used in references 22 and 33- 



(b) Relations used in present method. 

Figure 1.- Geometric relations associated with the application of the 

dovnvash boundary condition. 





Figure 2 . - Plan forms of wings for which flutter characteristics were calculated. 
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Figure 3.- Variation of flutter speed, with Mach number for wing For calculated points 

p = O.OO38OO slugs/cu ft. 
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Figure 5 .- Variation of flutter speed with Mach number for wing 44 5R. For calculated points 

p = 0.002378 slugs/cu ft. 






0.003500 slugs/cu ft. 
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Figure 8.- Variation of flutter speed with Mach number for wing b^>2. For calculated points 

p = 0.002700 slugs/cu ft. 
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Figure 9.- Variation of flutter speed with Mach number for wing 430. For calculated points 

p = 0.003700 slugs/cu ft. 
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Figure 11.- Variation of flutter speed with Mach number for wing 400R. For calculated points 

p = 0.003100 slugs /cu ft. 
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Figure 12.- Variation of flutter speed with Mach number for wing M+51- For calculated points 

p = 0.003200 slugs/cu ft. 
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Figure 14.- Variation of flutter speed with Mach number for wing 7001. 
For calculated points p = O.OO 55 OO slugs/cu ft. 
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Figure l6.- Variation of flutter frequency with Mach number for wing For calculated 

points p = 0.003000 slugs/cu ft. 




Figure 17.- Variation of flutter frequency with Mach number for wing 445R- For calculated 

points p = 0.002378 slugs/cu ft. 
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Figure 20.- Variation of flutter frequency with Mach number for wing kj2. For calc ula ted points 

p = 0.002700 slugs/cu ft. 
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Figure 22.- Variation of flutter frequency with Mach number for wing 1-00. For calculated points 

p = 0.002378 slugs/cu ft. 




Figure 23>- Variation of flutter frequency with Mach number for wing 1+00R. For calculated 

points p = 0.003100 slugs/cu ft. 
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Figure 24.- Variation of flutter frequency with Mach number for wing 4451* For calculated 

points p = 0.003200 slugs/cu ft. 
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Figure 25.- Variation of flutter frequency with Mach number for wing 4001. For calculated 

points p = 0.002378 slugs/cu ft. 




Figure 26 . - Variation of flutter frequency with Mach number for wing 7001. For calculated 

points p = 0.005500 slugs/cu ft. 
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(b) M =« 0.50. 
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(e) M = \J~2. 
Figure 28.- Concluded. 



(a) M = 0. 

Figure 29-- Distributions of static aerodynamic parameters for wing 645 . 
Symbols indicate -values of C^ a ^ n used in the flutter calculations. 
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(b) M = 1.15470. 



(c) M = \[ 2 . 


Figure 31* ~ Concluded. 
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(a) M = 0. 



(b) M = O.75. 

Figure 32.- Distributions of static aerodynamic parameters for wing 400. 

Symbols indicate values of C 7 used in the flutter calculations. 
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(d) M = l.JO. 
Figure Jk.- Concluded. 
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(a) M = 0. 
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(b) M = O. 75 . 

Figure 35*- Distributions of static aerodynamic parameters for 
wing 7001. Symbols indicate values of Cj used in the 


flutter calculations. 
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Figure 37*- First torsion and first and second bending mode shapes for 

a uniform cantilever beam. 
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Figure 41.- Circulation functions obtained from, the aerodynamic coeffi- 
cients given in reference 27 for two-dimensional airfoils oscillating 
in compressible flow. 
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